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ABSTRACT 


The  basic  requirement  for  the  development  of  a  six-degree- 
of-freedora  digital  simulation  model  are  outlined,  and  three  co¬ 
ordinate  systems  are  specified  which  are  both  adequate  and  con¬ 
venient  for  such  a  development.  The  nontrivial  coordinate  trans¬ 
formations  ar^  shown. 

The  equations  of  motion  are  developed,  with  indications  of 
the  standard  assumptions.  The  principal  forces  and  moments  are 
discussed.  The  thrust  misalignment  effects  are  derived  and  the 
formulas  for  the  jet  damping  moment  are  included. 

Aerodynamic  forces  and  moments  are  discussed,  and  the  necess¬ 
ary  aerodynamic  angles  are  defined.  These  forces  and  moments  are 
described  from  the  stability  derivative  point  of  view,  and  compon¬ 
ents  of  tne  major  forces  and  moments  are  derived. 
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INTRODUCTION 


Until  approximately  ten  years  ago,  the  study  of  rocket  tra¬ 
jectories  was  primarily  academic,  with  only  the  fundamental 
concepts  considered  and  very  liberal  assumptions  made.  Then, 
during  the  IGY  and  Tre-IGY  firings  of  unguided  soundirc  rockets 
it  was  discovered,  somewhat  embarrassingly  so,  that  simple 
"Kentucky  Windage"  'ith  the  carefully  moistened  finger  thrust 
aloft  was  not  adequate  to  determine  the  course  of  an  unguided 
rocket  and  that  a  more  objective  technique  was  required. 

Some  of  the  early  research  in  trajectory  simulation  and  wind 
effect  tabulations  for  unguided  rockets  was  performed  by  Lcwi.^ 
(1949),  Rachele  (1958),  and  Daw  (1958).  These  early  efforts, 
although  based  upon  restrictive  assumptions,  greatly  increased 
man’s  knowledge  about  the  problem  and  provided  techniques  for 
support  of  rocket  firings  which  were  used  for  several  years  with 
reasonable  success. 

The  ever- expanding  capability  of  the  electronic  computer 
has  augme  .ed  research  efforts  in  trajectory  simulations  by 
providing  the  computational  resources  required  for  precise  tra¬ 
jectory  simulations.  During  the  past  five  years  several  satis¬ 
factory  trajectory  simulation  models  have  been  developed;  some 
of  these  are  listed  in  the  references. 

There  have  been  two  unfortunate  occurrences  during  this 
rapid  advance  in  the  state  of  the  art:  (1)  in  the  majority  of 
the  literature  the  development  is  directed  toward  the  solution 
of  a  particular  problem  with  little  or  no  attention  given  to  the 
solution  of  the  basic  problem,  and  (2)  the  existing  literature 
has  been  published  as  technical  reports  of  limited  distribution 
which  are  difficult  to  obtain. 

This  paper  is  designed  to  outline  the  basic  considerations 
for  the  development  of  a  computer  program  for  trajectory  simula¬ 
tions.  The  discussion  is  not  pointed  toward  the  solution  of  a 
particular  problem;  however,  in  certain  instances,  the  techniques 
for  such  restrictions  are  outlined. 
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COORDINATE  SYSTEMS  AND  TRANSFORMATIONS 


The  number  and  nature  of  the  coordinate  systems  utilized  in 
a  simulation  model  depend  primarily  upon  the  complexity  of  simula¬ 
tion  requirements.  Three  coordinate  systems  are  presented  here 
which  are  adequate  for  most  trajectory  simulations. 

The  output  from  the  simulation,  as  well  as  the  initial  condi¬ 
tions,  is  usually  desired  with  reference  to  a  ground-fixed  launcher. 
This  system  will  be  referred  to  as  the  launcher  coordinate  system. 
The  only  convenient  system  in  which  to  compute  the  aerodynamic 
forces  and  moments  is  a  system  (commonly  called  a  body  system) 
which  is  affixed  to  and  moves  with  the  rocket. 


The  launcher  system  (denoted  X',  Y',  Z')  has  its  origin  at 
the  launcher  and  rotates  with  the  earth.  The  positive  X'  axis 
points  east;  the  positive  Y'  axis  points  north;  and  the  positive 
Z'  axis  points  'up'’.  The  body  system  (x,  y,  z)  has  its  origin 
at  the  center  of  gravity,  C  ,  of  the  missile.  The  x  axis  coincides 
with  the  longitudinal  axis  ind  is  positive  toward  the  nose.  Pre¬ 
cise  orientation  of  the  y  and  z  axes  is  somewhat  irrelevant  as 
long  as  the  system  remains  right-handed.  The  inertial  system  (X, 

Y,  Z)  has  its  origin  at  the  center  of  the  earth.  This  system  is 
oriented  so  that  the  X  and  Y  axes  lie  in  the  equatorial  plane  and 
the  Z  axis  is  coincident  with  the  earth's  axis  of  rotation  and 
positive  toward  the  North  Pole.  This  system  does  not  rotate  with 
the  earth.  Although  the  exact  orientation  of  the  X  and  Y  axes  is 
somewhat  arbitrary,  it  is  convenient  to  define  one  of  these  such 
that  it  initially  passes  through  either  the  longitude  of  the  launch¬ 
er  or  through  longitude  0. 


A  linear  transformation  between  any  two  of  these  systems  will 
be  denoted  by  T„_  where  the  '  ft -hand  subscript  denotes  the  domain 
of  the  mapping. x^§ince  the  coordinate  systems  are  all  orthogonal, 
the  inverse  of  the  transformation  is  just  the  transpose.  The  trans¬ 
formation  TX21,  is  easily  determined  by  geometric  and  trigonometric 
consideration!.  This  transformation  depends  upon  the  earth-model 
considered  (e.g.  spherical,  oblate  spheroid,  pear-shaped,  etc.)  and 
the  earth's  rotation. 


The  ransformation  T„2  is  not  so  easily  obtained.  Two  methods 
will  be  given  below  for  ooctining  this  transformation.  The  develop¬ 
ment  of  these  methods  is  lengthy  and  will  be  omitted.  The  first 
method  is  based  on  Euler  angles,  and  a  development  can  be  found  in 
much  of  the  literature,  e.g.,  Lass  (19S0).  The  second  method  is  based 
on  direction  cosines;  a  development  has  been  presented  by  Duncan  (1966) 
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The  x-y  plane  will  intersect  the  X'-Y*  plane  in  a  line, 
called  the  nodal  line  N.  Let  Q  be  the  angle  between  the  z  and 
Z'  axes,  ip  the  angle  between  the  X*  and  N  axes,  and  4>  the  angle 
between  the  N  and  x  axes.  (This  is  but  one  of  several  ways  by 
which  the  Euler  angles  can  be  defined.)  The  angles,  ip,  0,  $ 
completely  specify  the  relative  orientation  of  the  two  systems; 
hence,  the  rotation  matrix,  the  matrix  of  T^2x,  can  eas^^Y  be 
determined  once  the  values  of  these  angles  are  determined. 

Let  w  define  the  angular  velocity  of  the  x,  y,  z  system 
relative  to  the  X,  Y,  Z  system.  Then 


w  »  pi  ♦  qj  ♦  rk, 


(1) 


where  i,  j,  and  k  are  unit  vectors  along  the  x,  y,  z  axes.  It  is 
shown  by  Lass  (1950)  that. 
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If  the  Euler  angle  technique  is  used  to  determine  the  trans¬ 
formation  T„-  »  the  system  of  equations  (2)  becomes  three  of  the 
equations  or  motion  in  the  simulation  model. 

The  second  technique  employs  differential  equations  involving 
dt  and  the  direction  cosines.  Let  (£. ,  ^2>  £  ),  (m  ,  a„,  m3)  and 
(n.,  n_,  n.)  be  the  respective  direction  cosines  or  thfc  x,  y,  and 
z  axes^in  the  X,  Y,  Z  system.  Then  the  matrix  of  TX2x  is 


n  i  r 

■  1  2  3 


B1  m2  b3 


t 

t 


n2  n3 


(3) 
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It  has  been  shown  by  Duncan  (1966)  that 


m. 

1 


n. 

1 


rm.  -  qn.  i  *  1, 

l  n  l 

pnt  -  rf.  i  «  1, 

qt  -  pm.  i  “  1  > 


2,  3 
2,  3 
2,  3 


(4) 


If  the  second  method  is  used  to  determine  the  transformation 
T„2  ,  then  the  system  of  equations  (4)  becomes  nine  of  the  equations 
or  motion  in  the  simulation  model. 

Both  of  the  methods  discussed  above  have  certain  undesirable 
features.  In  method  one  it  is  possible  for  one  of  the  Euler  angles 
to  become  undefined;  method  two  requires  the  integration  of  a  larger 
system  of  equations. 


THE  EQUATIONS  OF  MOTION 


The  equations  involving  the  moments  of  inertia,  aerodynamic 
forces  and  moments,  and  the  thrust  forces  are  greatly  simplified 
if  expressed  in  the  body  coordinate  system;  hence,  this  system  will 
be  used.  The  two  basic  equations  which  define  the  motion  of  a  rigid 
oody  are 

F  ■  (mV) ,  (5) 

and 

Mr  *  (i5)  (6j 

where  F  represents  the  external  forces, 

represents  the  external  moments, 

H  represents  the  angular  momentum. 


Numerical  analysis  of  these  vector  equations  requires  their 
resolution  into  vector  components  and  definition  of  the  scalar 
coefficients.  These  manipulations  are  discussed  in  detail  in  many 
texts  in  Mechanics.  The  essential  steps  are  reviewed  here,  however, 
for  completeness. 

To  determine  the  translational  acceleration^  consider  a  point 
P  defined  in  the  (x,  y,  z)  system  by  the  vector  r.  Let  the  origin 
of  the  (x,  y,  z)  system  with  respect  to  the  X,  Y,  Z  system  be 
specified  by  the  vector  R. 

Then 


r  •  xkx  ♦  yky  ♦  zkz, 

$  «  Xkx  ♦  Yky  ♦  Zkz.  (7) 


In  the  development  of  the  equations  of  motion  the  forces 
are  usually  assumed  to  act  through  the  center  of  gravity  of  the 
rocket,  i.e.  r  «  0.  It  is  more  convenient  to  express  the  velocity 
of  the  origin  of  the  body  system  in  body  axis  coordinates;  these 
are  usually  denoted  by  u,  v,  w. 

The  equation  F  •  is  true  in  the  X,  Y,  Z  system  but  must 
be  modified  if  the  forces  are  to  be  computed  in  the  body  system. 

To  this  end,  let  S  be  a  vector  in  the  (x,  y,  z)  system.  Then 


DS  dS 

— -  ■  -! -  ♦  U 

dt  dt 


S 


(8) 


where  the  symbol  -rr  indicates  diff jrentiation  in  the  (X,  Y,  Z) 
system  and  ^  refers  to  the  body  system. 

Now  the  first  equation  of  motion  becomes,  in  component  form: 


F 

F 

F 


x 

y 

z 


m(u  -  rv  ♦  qw) , 
m(v  -  wp  ♦  ur), 
m(w  -  uq  ♦  vp) . 


(9) 


5 


These  equations  determine  the  translational  motion. 


The  equations  expressing  the  rotational  motion  are  obtained 
in  a  straightforward  manner.  The  angular  momentum  of  a  body 
about  its  center  of  gravity  is  given  by 


H  «  [I  p  -  I  q  -  I  r]k  ♦  [-1  p  ♦  I  q  -I  r]k 
u  xxr  xyM  xz  x  u  xy^  yy^  yz  y 

♦  [-1  p  -  I  q  +  I  r]k  .  (10) 

*•  xzy  y?M  zz  J  z  v  * 


It  is  the  general  practice  at  this  point  in  the  derivation 
of  the  equations  of  motion  to  assume  that  the  body  axes  are 
principal  axes  of  inertia.  Under  this  assumption  one  has 


it  »  I  pk  +  I  qk  +  I  rk  . 
xx*  x  yyM  y  zz  z 


(ID 


Hence, 


H  *  (I  p  +  I  p)k  +  [i  q  +  I  q]k  +  [i  r  +  i  r]k 
v  xx v  xxr  x  yyn  yyMJ  y  zz  zz  z 


♦'Ik+Ik+Ik 
xx  x  yy  y  zz  z 


(12) 


Now  k  *wxk  ,  k  ■  w*k  ,  k  =  u>xk  .  Hence 
x  x’  y  y  ’  z  z 


H  -  [I  p  +  I  p  +  (I  -  I  )qr]k 
4  xxr  xxr  v  zz  yy  x 


+(l  q+I  q  ♦  (I  -I  )pr]k 
yy  yyn  xx  J  y 


(13) 


+[l  r  ♦  I  r  +  (I  -I  )pq]k  . 
*■  zz  zz  v  yy  xx  z 


The  x,  y,  z  components  of  the  total  external  moment  are  commonly 
known  as  L,  M,  N,  respectively.  Now  the  rotational  equation  of 
motion  becomes,  in  component  form. 
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L 


*  I  p  ♦  I  p  +  (I  -  I  )qr 

xxr  xxr  zz  yy'^ 

•  • 

M  =  I  q+I  q  ♦  (I  -I  )pr 

yy^  yy '  xx  zz  r 

N  *  I  r  ♦  I  r  ♦  (I  -  I  )pq, 

zz  zz  yy  xx/tn 


(14) 


FORCES  AND  MOMENTS 


The  forces  and  moments  acting  on  a  rocket  are  due  to  three 
specific  effects  These  are  the  thrust,  the  gravitational  attrac¬ 
tion,  and  the  aerodynamic  features  of  the  rocket.  The  forces  and 
moments  are  expressed  in  the  body  coordinate  system. 

Gravitational  Effects: 

The  gravitational  force,  mg,  is  $asy  to  compute.  However, 
the  exact  magnitude  and  direction  of  g  depend  upon  the  earth  model 
chosen  for  the  simulation.  Since  this  subject  is  discussed  quite 
thoroughly  in  the  literature  it  will  not  be  discussed  here. 

Thrust  Effects; 

Let  m  be  the  mass  of  the  rocket  including  the  unspent  fuel 
and  let  Am  be  the  change  in  mass  (due  to  burning  of  fuel)  during 
a  small  time  interval  At.  By  the  law  of  conservation  of  momentum, 
the  momentum  at  time  t  is  equal  to  that  at  time  t  ♦  At. 


m\/  =  (m  +  Am)  $  *  A\f)  ♦  Am(\^  -  \f) 


(15) 


where  V  is  the  velocity,  in  the  body  system,  of  the  exit  gases. 
Hence,  ^  *  -  m$  .  This  is  the  force  on  the  rocket  due  to  the 
changing  momentum. 

Besides  m^e  there  is  an  additional  force  due  to  the  differ¬ 
ence  between  the  pressure  at  the  exit  nozzle  and  the  atmospherifc 
pressure.  If  the  pressure  at  the  exit  nozzle  is  P  ,  the  atmos¬ 
pheric  pressure  P  ,  and  the  area  of  the  exit  nozzll  A  ,  then  this 
additional  force  Is  A  (P  -  P  ) ,  giving  a  total  thrusf  of 

0C  8. 
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(16) 


3b 


t 


"v.  ♦  -  W 


If  the  thrust  is  measured  at  a  test  stand  at  an  atmospheric 
pressure  P  it  would  be 


♦  A  (P 
e  e 


Hence, 


T  *  ♦  A  (P  „  -  P  ) 

s.t.  e  s.t.  a' 


(17) 


(18) 


Since  a  rocket  rotates  about  a  transverse  axis  during  burning, 
the  gases  must  be  accelerated  laterally  as  they  flow  through  the 
nozzle.  This  lateral  acceleration  produces  the  so-called  jet  damp¬ 
ing  moment.  The  following  expression  for  the  jet  damping  moment 
was  derived  by  Brown  et  al.  (1961): 


M  *  m(f?  -  f?)w  , 
P  J  P  P 


(19) 


where  m  is  the  mass  flow  rate,  w  is  the  instantaneous  pitching 
velocity,  l.  is  the  distance  between  the  vehicle’s  C  and  the  exit 
nozzle,  and^-f  is  the  distance  between  the  C  of  the8vehicle  and  the 
propellant  8 

The  components  of  this  moment  are,  for  a  symmetric  rocket. 


M.  .  m{l2  -  l2)q, 
J  J  P  ^ 


**.  «  m(4?  -  ^)r. 
J  J  P 


(20) 


The  rocket  thrust  is  capable  of  producing  components  of  force 
and  moment  along  each  of  the  body  axes.  These  may  be  due  to  a  mis¬ 
alignment  of  the  thrust  vector  with  respect  to  the  x-axis  or  to  an 
off-center  installation  of  the  rocket  motor.  Let  x.,  y,,  and  z^  be 
the  Xj  y,  z  coordinates  of  the  i-th  exit  nozzle.  Lit  T^  be  the 
thrust  vector  of  the  i-th  exit  nozzle.  Suppose  the  i-th  thrust  vector 


8 


is  oriented  as  shown  in  Fig.  1.  Then  the  components  of  the  i-th 
thrust  vector  are 


T.  =  T.  cosA_ 
ix  1  T. 

1 


T.  «  T.  sinA_  cos^ 
ly  l  T.  *T. 


T  -  T  sinX  sin*T 
1  1 


(21) 


The  components  of  the  total  thrust  vector  f  *  ][  f. 

i 

are  T  *  J  T.  ,  T  *  T  T.  ,  T  *  J  T.  .  It  follows  immediately 

x  v  ?x  y  v  iy’  z  4-  iz  7 

1  J  1  J  1 

-+•  -* 

from  the  relation  M.  =  r.  x  F.  that  the  components  of  the  i-th 


moment  are 


U  =  y.T.  -  z.T. 

i.  ri  1Z  11 
l 


M_,  =  z.T.  -  x.T 
T.  1  IX  1 


1Z 


(22) 


N_  =  x.T.  -  y.T. 
T^  l  ly  7i  ix 


where  r^  *  (x^,,  y^,  z^)  and  ,  and  NT  are  the  i-th 

thrust  moments  about  the  x,  y,*and  i  axes,  respectively.  The 
components  of  the  total  thrust  moment  M  *  £  NL  are 

LT“lLT'MT-K.-NT*?NT.'  <23> 

i  i  ii  ii 


Aerodynamic  Forces  and  Moments 

There  are  several  angles  which  are  used  to  calculate  the 
aerodynamic  forces  and  moments.  These  angles  are  shown  in  Figure 
2.  They  can  be  expressed  in  terms  of  the  velocity  components  as 
follows.  Let  w  ,  w  ,  w  be  the  x,  y,  z,  components  of  the  wind. 
The  components  8f  tKe  velocity  of  the  rocket  relative  to  the  wind 
are  u*  «  u  -  w^,  v'  ■  v  -  w^,  w*  ■  w  -  w^  and  the  relative  speed 

is  V  *  [(u')2  +  (v ' ) 2  ♦  (w*)2]5*.  The  angle  of  attack,  a,  the 

a 
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angle  of  sideslip,  8,  and  the  absolute  angle  of  attack,  5,  are 
defined  by 


.  -IrWN 

tan  L — rJ 


tan"1  [ 


3 


(24) 


/(u*)2  ♦  (w')2 


tan'1  ’)2] 


The  auxiliary  angle  of  attack,  a*,  and  the  auxiliary  angle  of 
sideslip,  8*,  are  given  by 


*  “1  r 

a  =  tan  L 


w' 


An')2  +  (v')2 


3 


(25) 


*  -I  rV'i 

8  *  tan  [—] 


The  general  technique  for  specifying  the  aerodynamic  forces 
and  moments  utilizes  the  concept  of  stability  derivations.  Stability 
derivations  have  been  discussed  in  considerable  detail  by  Neilsen 
(1960).  The  basic  ideas  and  the  application  of  the  concept  will  be 
presented  below. 

The  formulas  for  computing  the  forces  and  moments  are 


F  »  C^'S 

Mj.  ■  C^q’Sd  (26) 


where  and  C.  are  dimensionless  coefficients,  q'  is  the  dynamic 
pressure,  S  ana  d  are  the  reference  arsa  and  reference  length, 
respectively.  The  above  equations  are  usually  written  in  compon¬ 
ent  form.  They  become 
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Fx  •  Cx«’S’ 
Fy  ’  Cyq’S, 

F,  =  C2q'S, 

M  -  C  q’Sd, 
m 


N  -  Cnq'Sd, 


(27) 


L  -  C^q'Sd, 


where  C  ,  C  ,  C  ,  C  ,  C  ,  and  C i  are  considered  to  be  functions  of 
severalxvar¥ablls.  ^testability  derivations  are  simply  partial 
derivatives  of  these  functions.  The  procedure  by  which  stability 
derivatives  are  applied  is  best  described  by  an 


example.  Suppose  *  ffc^,  a2»  •••»  an)  •  Let  cmo  *  af/3a*  an<* 

suppose  C  is  known  for  some  value  (a,  .  ...,  a  );  call  this 

value  C  t  Then  if  each  of  the  a.'s  dnanges  byna  small  amount  da. 
mo  3  j 

n 


C 


m 


(28) 


It  is  usually  assumed  that  (a.  ,  a.  ,  ...,  a  )  ■  (0,  ...  0)  and 
that  C  is  linear  for  all  realizable  neighborhoods  of  this  point. 
Under  Wlis  assumption  da^  is  approximated  by  a^  and  the  above  equa¬ 
tion  becomes 


n 

C  «  C  ♦  I  a.C 
m  mo  1  ma^ 


(29) 


The  Forces  and  Moments  due  to  Air  Resistance 


The  specific  stability  derivatives  included  in  the  development 
of  a  simulation  model  depend  upon  the  purpose  of  the  computation  and 
the  details  required  therein.  (The  availability  of  numerical  values 
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for  the  stability  derivatives  is  often  another  controlling  factor.) 
Since  a  derivation  of  the  applications  of  all  the  various  combina¬ 
tions  of  the  stability  derivatives  defined  by  Neilsen  (1960)  would 
lead  to  voluminous  formulations  of  questionable  value,  only  those 
resulting  from  the  effect  of  air  resistance  will  be  discussed  here. 

The  force  due  to  air  resistance  is  broken  into  components 
parallel  to  and  perpendicular  to  the  x-axis;  these  components 
are  referred  to  as  the  axial  and  normal  force,  respectively.  The 
standard  notations  for  the  dimensionless  coefficients  are  C.  and 
Cj.,  (The  force  is  sometimes  resolved  into  a  different  reference 
frame  and  the  components  are  referred  to  as  drag  and  lift.)  The 
normal  force  lies  in  the  plane  of  the  x-axis  and  the  V  vector; 
the  direction  is  such  that  the  force  tends  to  decreased. 

Most  rockets  possess  a  property  called  90-degree  roll  symmetry. 
This  means  that  the  physical  characteristics  of  the  rocket  remain 
unchanged  if  the  vehicle  is  rotated  90°  about  the  x-axis.  Thus  for 
the  symmetric  vehicle,  C  *  C  and  C  «  C  ;  hence,  the  stability 
derivatives  of  these  functions  are  eSfual.n 

The  coefficient  C  is  primarily  a  function  of  the  angle  6.  If 
it  is  assumed,  and  it  Sften  is,  that  C  is  a  function  of  6  alone 
then  it  is  easy  to  see  from  Figure  2  tKat  the  components  of  the 
normal  force  and  the  moment  contributions  due  to  this  force  are 


B, 

«\  (30) 


o  . 


6 


where  C  i»  the  moment  coefficient  which  results  from  the  force. 

M 

The  rotation  of  the  vehicle  generates  a  damping  moment  (due 
to  the  air  resistance).  The  standard  formulation  of  the  stability 
derivatives  is: 


-  C  sin 
na 


-  C  sin 
na 


♦  C  sin 
not 


c_  • 


-  sin 
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(31) 


%  *  ,c«'*(Sr) 

c«q  •  iC«/s(§r) 


c«  ’  *VJ(Kr> 


The  above  formulations  give  the  following  expressions  for 
the  forces  and  moments  due  to  air  resistance  under  the  assump¬ 
tion  of  a  90°  roll  symmetric  vehicle: 


*  c,«'s 


"y  *  "Cna  5ln 


-C  sin  a  q'S 
na 


(32) 


L 

H 

N 


[C^  ♦  Ct  (J$-)]q’S  d 
[C^  sin  a*  ♦  C^  (^-)]q'S  d 


[-C^  sin  8  ♦  CBr  (j$-)]q'S  d. 


CONCLUSIONS 


The  principal  considerations  in  the  development  of  a  six-degree- 
of- freedom  trajectory  simulation  model  have  been  discussed.  Although 
any  particular  simulation  model  may  include  several  factors  not 
discussed  in  this  paper,  it  must  include  ai  1  follow  the  basic  prin¬ 
ciples  outlined  I.erein.  This  treatise  should  provide  a  suitable  back¬ 
ground  to  the  researcher  who  is  required  to  develop  a  simulation  model 
and  should  indicate  how  the  various  idiosyncracies  of  his  particular 
problem  may  be  handled. 
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